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Abstract. - We study a model of an active gel of cross-linked semiflexible filaments with additional 
active linkers such as myosin II clusters. We show that the coupling of the elasticity of the 
semiflexible filaments to the mechanical properties of the motors leads to contractile behavior of 
the gel, in qualitative agreement with experimental observations. The motors, however, soften the 
zero frequency elastic constant of the gel. When the collective motor dynamics is incorporated in 
the model, a stiffening of the network at high frequencies is obtained. The frequency controlling the 
crossover between low and high frequency network elasticity is estimated in terms of microscopic 
properties of motors and filaments, and can be as low as 10 _3 Hz. 
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Introduction. — The mechanical properties of cells 
control many biological functions, including the sensing 
and generation of forces, cell motility and cell division. 
The response of the cell to mechanical stimuli is mediated 
by the cytoskeleton, a network of semiflexible filaments 
(F-actin, microtubules and intermediate filaments) linked 
by a variety of passive and active proteins. [1,2] The cy- 
toskeleton is maintained out of equilibrium by chemical 
reactions that drive force generation by motor proteins, as 
well as by filament treadmilling. A variety of recent exper- 
iments have measured the remarkable rheological proper- 
ties of this intrinsically nonequilibrium polymer network. 
These include bulk and microrhcology of in vitro stabilized 
networks of cytoskclctal filaments with a controlled con- 
centration of various crosslinkers, as well as in vivo whole 
cell rheology. 

Cross-linked entangled actin networks are viscoclas- 
tic solids, with a time-dependent mechanical response 
(stress a) to deformation (strain 7). These networks have 
both viscous and elastic responses characterized by loss 
G"(u>) ~ cr/7 and storage moduli G'(oj) ~ cr/7, respec- 
tively. For cross-linked gels, the elastic (storage) modu- 
lus dominates the mechanical response and reaches a fre- 
quency independent plateau Go at low frequencies (less 
than 1Hz). Experimentally Go is found to depend strongly 



on cross-link density and can vary from 0.1 - 100 Pa [3]. 
For frequencies above 1Hz, both the storage and loss mod- 
uli show a high frequency behavior G', G" ~ w 3 / 4 charac- 
teristic of semiflexible polymer dynamics [4] . 

Measurements of the mechanical properties of cells 
yield, however, quite different behaviour [7]. The low 
frequency (< 10Hz) shear moduli are observed to be- 
have as, G',G" ~ G*(u>/ui*) a , with a small exponent 
a ~ 0.15 - 0.2, G* ~ 10 2 - 10 3 Pa and uj* ~ 1 Hz [8-13]. 
Significantly, the magnitude of G* is much higher than 
the typical plateau moduli of purified in-vitro actin gels. 
While increasing cross-linker density can significantly en- 
hance the elastic modulus [3] , it is surprising that it would 
have such a dramatic effect on the loss modulus. It was 
recently suggested that the remarkable stiffening of the 
low frequency linear response of active gels may be due 
to the internal stresses generated by the presence of ac- 
tive crosslinkers, such as myosin II minifilaments [14,15]. 
Recent quantitative experiments studying the mechanics 
of in-vitro networks of F-actin, with passive (a— actinin) 
and active (muscle myosin II) cross-linkers, have shown 
both stiffening [14] and contractile behaviour [17] of these 
reconstituted networks. Interestingly the contractile be- 
haviour has been shown to appear only in a narrow con- 
centration range of passive cross- linkers. 
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In this letter we present a theoretical description of ac- 
tive gels which can explain both the contractile behaviour 
and the intermediate- frequency stiffening of these systems. 
The minimal element from which the active gel is con- 
structed is a pair of filaments cross-linked by an active 
cluster of molecular motors. This has been a useful start- 
ing point for explaining the properties of soft active mate- 
rials in both the fluid [18] and the gel phase [15]. Our work 
shows that the coupling of the elasticity of the semiflexiblc 
filaments to the motor dynamics plays a crucial role in con- 
trolling the rheology of the network. We find that active 
clusters lead to contraction of the gel which has a more 
dense ground state than a gel with the active crosslinks 
replaced by passive ones. The zero frequency stiffness of 
our model active gel, when perturbed from this ground 
state, is lower than that of the corresponding passive gel, 
in apparent contradiction with experiments. However, at 
higher frequencies the collective dynamics of the motors 
stiffens the gel as compared to the passive case. This qual- 
ititivc behaviour is obtained both in the regime of linear 
chain elasticity and when taking account of nonlincaritics. 

Model. — We consider an ideal semiflexible polymer 
network with both permanent and active crosslinkers. The 
network consists of isotropically oriented stiff polymer seg- 
ments of length shorter than their persistence length sub- 
jects to rigid constraints due to the permanent crosslinks. 
The motor clusters act as dynamic cross-linkers which ap- 
ply equal and opposite forces to pairs of filaments. We do 
not consider the effect of entanglements. 

We parametrize each filament by a curve R(s), with 
< s < L and L the contour length. We consider small de- 
viations from a straight configuration of the polymer seg- 
ment and decompose deformations of length scales smaller 
than the persistence length L p = n/ksT, with k the bend- 
ing rigidity, in transverse and longitudinal components by 
writing R(s) = i?||(s)u + r±(s) = [s — r\\(s)]u + r±(s), 
where u is a unit vector giving the orientation of the seg- 
ment and r ■ u — [16]. In a cross-linked gel, the free 
energy of each filament is given by 



F = 



His; 



R| 2 -o-ou-d s R\ 



(1) 



The first term is the usual bending energy of a worm-like 
chain and the second one describes the tension <jq 0I the 
filaments due the permanent crosslinkers. The filament 
satisfies the boundary conditions R(0) = 0, R(£) = L n u 
and [<9 2 R(s)] s=0 = [<9 2 R(s)] s=L = 0, where L is the dis- 
tance between fixed (passive) crosslinks [19]. Transverse 
and longitudinal deviations are coupled by the constraint 
of inextensibility, d s r\\ = ±|<9 s rJ 2 + C(|9 s r ± | 4 ) [16]. 

The effective longitudinal response function of a fila- 
ment is evaluated by averaging over the transverse fluctu- 
ations, with the result 



<a s r||) Q = — { XC ° J } = F(L,o- ,k), 





Fig. 1: Top: A schematic representation of two semiflexible 
filaments crosslinked by an active myosin cluster (shown in 
red) and linked to neighboring filaments (blue dashed lines) by 
passive linkers (blue dots). Bottom: The minimal elastic ele- 
ment of our active gel, consisting of two antiparallel filaments 
crosslinked by a motor cluster. The motor cluster is modeled 
as a spring of stiffness k m that exerts equal and opposite forces 
of magnitude / on the two filaments. The pair of filaments is 
maintained at a fixed distance Lq by passive crosslinkers. The 
asymmetric shape of each filament indicates that due to their 
nonlinear elastic properties the filaments are easier to compress 
than to stretch. 



and x = LJ 1 ^-. The end-to-end length of the filament is 



L- 



d»<fl.r u > = L-[(r||(L)> -<ru(0))J 



(3) 
In the limit uq — ► the filament is roughened by thermal 



fluctuations and T ~ 
vcrsely, for x ^$> I T 



(2) 



^P [1 - x 2 /15] for x < 1. Con- 
~ fcsT/2^/KfTo and all wrinkles are 
pulled out by the applied tension so that L = Lq. 

Now we consider the effect of additional active cross- 
links on the mechanical properties of the gel. If the gel is 
kept under constant external tension, the activity of the 
motors changes the end-to-end distance L . Conversely, 
if Lq is changed by applying an external deformation, the 
active cross-links induce an additional tension on the fil- 
aments. Both types of response may be studied by eval- 
uating the change in extension of the filaments upon in- 
creasing of the local tension from o~o to erg + ct(s). The ad- 
ditional tension may be thought of as arising either from 
motor activity or from externally applied forces. To de- 
scribe this response we define an 'elastic' deformation field 
u ( s ) = [(-R||( s )) — {-R|| (s))o] , where (i?||(s)) denotes the 
effective longitudinal response to the total tension <7o + a . 
In the limit of large ctq, this yields a general relationship 
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between the deformation and the tension, 
d s u{s) = G[a ,<j(s)] , 



force balance equations 
(4) Co<9 t si = -(od t s 2 = -k m {si -s 2 + ui -u 2 ) -/, (6) 



with G = F{L,<7q,k) - P{L,a + a{s),K). Eq. (4) de- 
scribes the nonlinear elasticity of a semifiexible filament 
under tension and is the starting point of our analysis. 

The top image of Fig. 1 shows a schematic of a motor 
cluster crosslinking two semifiexible filaments, which are 
in turn bound at their ends by permanent crosslinks. The 
motors in the cluster walk towards the plus end of each fil- 
ament indicated by the double arrows, exerting equal and 
opposite forces ±/ on the two filaments [14,15,18], result- 
ing in additional tension. Such a crosslinkcd filament pair 
is the fundamental elastic unit in our model of an active 
gel. For an isotropic gel, with uncorrelated orientations 
of filament pairs and motor clusters, the elastic properties 
of the gel may be obtained by suitable angular averages. 
All the essential physics can, however, be obtained from 
the simplified one-dimensional model depicted in the bot- 
tom part of Fig. 1. Taking account of the orientation of 
the filaments only changes the numerical prefactors. Also 
we consider only anti-parallel filaments as motors do not 
generate significant tensions on pairs of parallel filaments. 
Our goal is to calculate how the end-to-end distance and 
the mechanical response of the filament pair are changed 
by the presence of the active crosslink. The force balance 
equations for the two filaments are 



of - erf + f = , a? - o\ - f = 



(5) 



R.L 



where er"' 2 ", for i = 1,2 are the additional tensions (on 
top of the "bare" tension ao) at the ends of each filament 
due to active cross-links and/or external forces (Fig. 1). 

Motor clusters. — We consider the response of the 
system on time scales long compared to the relaxation 
time of the longitudinal modes of the filaments so that 
we can ignore the dynamics of the deformation, w(s). We 
include, however, the finite-frequency response of the mo- 
tors, which is important when considering the response of 
the gel to frequency-dependent deformations. To do this 
we use the model introduced by Jiilicher and Prost [20], 
where a motor cluster actively sliding along a polar fila- 
ment is described as a collection of N motors rigidly at- 
tached to a backbone and moving along a polar periodic 
track. Under the action of both thermal and ATP-driven 
excitations, each motor in the cluster undergoes transi- 
tions between a strongly bound and a weakly bound state. 
We model our active cross-link as two such motor clusters 
linked by a spring of stiffness k m that couples the mo- 
tion along the top and bottom filaments. We denote by 
Si the position of the cluster along the z-th filament and 
by Ui = u(si) the elastic deformation of the i-th filament 
at that point and assume that the center of mass position 
■ s o = ( s i + s 2 )/2 of the motor cluster does not move. 

The dynamics of the points of attachment of the motor 
cluster interacting with the filaments is described by the 



where £o is a friction. The first term on the right hand 
side of Eq. (6) represents elastic forces within the motor 
cluster. In the last term, — / is the force exerted by the 
motors on filament 1. In a steady state, where the velocity 
of the motor cluster vanishes, / equals the stall force f s . 
As in Ref. [20], if the motor clusters are not exactly at 
stall force we expand the force up to first order in the mo- 
tors velocity and introduce an active frequency-dependent 
friction C,actiy) on the motors, which can be negative. For 
a sake of simplicity, we will assume here that the motor 
clusters have no spontaneous oscillations, which is the case 
if the ATP concentration is low enough. 

Static response. — We first examine the static re- 
sponse of a cross-linked filament pair to a change in the 
end-to-end distance Lq. We apply a force 2F to the fil- 
ament pair and calculate the resultant extension of the 
right-hand side of both filaments A(F). A value A(0) < 
at zero applied force corresponds to contractile behaviour. 

The conformation of each filament is described by a 
displacement Ui(s), with d s Ui = G [ao, o"i(s)]. These 
equations must be solved with with boundary conditions 
Ui(0) = and Ui(L) = A, for i = 1,2. Force balance on 
the filaments indicate that the tension has a jump discon- 
tinuity at the point of motor attachment. This implies 
that the derivative d s Ui(s) is also piecewise constant, with 
a jump discontinuity of magnitude / at s^. Requiring the 
displacement of each filament Ui = Ui(si) at the point of 
attachment to be continuous, we obtain 



A -Ui 
L - Si 



G[a ,a + a t ] , 

: G[(T0,(T0 + 0~i ] , 



(7) 
(8) 



for i = l,2. Eliminating u\ and u 2 from Eqs. (7) and (8), 
and from the stall condition, f s = f = — k m (As + u\ —U2), 
we obtain 



siG(of ) + (£ - si)G((r?) = A : 

s 2 G(a 2 L ) + (L-s 2 )G(a«) = A. 

( Sl - s 2 ) + Sl G{a^) - s 2 G(a%) 



h ri 



(9) 
(10) 

(11) 



where G(a) = G[<to,<tq + a]. From the two force bal- 
ance equations, Eqs. (5), it is evident that only two of the 
four tensions a, t ' are independent. It is convenient to 
eliminate two of the unknowns by introducing new forces 
^i — I ( (T 4 fl + <7 f): so that the force balance equations 
are automatically satisfied. Letting si ]2 = s$ ± As/2, 
F = (Ft + F 2 )/2 and 6 F = F ± - F 2 , Eqs. (9-11) yield a 
set of three coupled equations in three unknown As, SF 
and either A of F. The equations can be solved to obtain 
either the displacement A as a function of the total force, 
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A(F), or F(A). The solution will depend parametrically 
on the center of mass position of the motor cluster, sq. 

We solve Eqs. (9), (10) and (11) taking into account the 
nonlinear elasticity of the filaments, in the perturbative 
limit where all the motor induced forces are small com- 
pared to the bare tension <7o and solve the equations per- 

D r 

turbativcly in /. If all tensions a i ' are small compared 
to (To one can approximate G[<7o,<t] = ^2 n G n a n /nl ss 
Go + G"a 2 /2 + G"V 3 /6 + . . ., where G n = G n (er ) = 

-G[a ,a +a]\ ^^ Q , = ( d G[(J Q , o] / da) ^ CtC . 



da" 



(7=0 



To linear order in the total force F, we write 
A(F) = A(0) + -!— + 0(F 2 ) . 
The ground state deformation is given by 



A(0) 
L 



2 



G' 
Km Li 



;i-0)(G"-2(G') 2 ) 



(12) 



(13) 



where 4> = s /L. Since G'(cr ) > and G"(a ) < 
the ground state deformation is always negative, corre- 
sponding to a contractile system. This result is easily 
understood if we consider the limit where G" = for 
n > 1 and the filaments behave as a linear springs of 
elastic constant fco = (LG') _1 . In this case the ground 
state deformation can be written as A = /G'Aso/2, where 
Aso = — |t — 2(f>(l — (f>)-j£- is the ground state value of sep- 
aration of the motor clusters between the two filaments. 
Then A is easily obtained by equating the change in elastic 
energy when the filaments are stretched from L to L + A, 
given by 2[|fc (i + A) 2 - ±fc £ 2 ] ~ 2k Q LA to the work 
fAs done by the motor clusters on the filaments. The 
effective stiffness of the network is given by 



1 

k cff 



1 



l\3 



+ Lf4>{\ - cf>) G'"/2 - 2G'G" + (G) 



.JL 



\G" + (G') 



/\21 



(14) 



The average ground state deformation and the effective 
stiffness of the element are shown in Fig. 2 as functions 
of L, with A(F) = L dcj)A(F). Contractile behaviour is 
observed for all L and vanishes as L — > 0, reflecting the 
higher resistance of short filaments to compression. On the 
other hand, the active crosslinks always decrease the zero 
frequency stiffness of the gel which vanishes as L — > oo. 

Our perturbative analysis captures the qualitative ex- 
perimental observation of contractile behaviour, but yields 
softening of the gel in contradiction with experiments. To 
address this we incorporate in the next section the finite 
frequency response of the motor clusters. 

Finite frequency behaviour. — We now consider 
the finite frequency behaviour near stall and take into ac- 
count the collective dynamics of the motors within an ac- 
tive crosslink cluster. We consider time-scales long com- 
pared to the relaxation time of the filaments, but include 
the finite- frequency response of the motors. 




Ktt 


/k 














|.2 4 


6 


8 


10 


0.98 












0.96 
0.94 




V 








0.92 

0.9 

0.88 












0.86 




\ \ 









f Wao Ik 



Yq Wao Ik 



Fig. 2: Plot of — A(0)// 2 !/ (top) and zero frequency stiffness 
fc e g- (bottom) for l/fc m = 0, ■£- = 0.4 as a function of L^/oo/k 

for 3 values of L p y/a /K =0.2(red), 0.17(blue), 0.13(green). 



We apply an oscillatory deformation to the end-points 
of the right side of the filaments, A = Ao + eAi(i), while 
keeping the left side fixed. Here Ao = A(0) is the static 
displacement at zero external force, corresponding to the 
motors' stall force. This perturbation will result in motors 
sliding along the filaments. To discuss the response of 
the system we work in Fourier space and linearize in all 
the deviations /—>■/ + Sf^ , As — ► As + 5s u , F — > 
F u , 8F -> SF + SF^ ,where /, As, F and SF denote 
static quantities at stall condition as defined earlier. From 
equation (6) the dynamics of the relative displacements of 
the two motor cluster's heads is given by 



iiuCoSs u 



ZK r 



1 



l + TT[GK L )+G(a 2 L )] \Sa u -25f l 



-2k m [ Sl 8otG'(^) - s 2 5atG'(a L )} (15) 



where <5< = F^ + ^SF^ + SL) and Sa^ = F UJ -±{5F U + 
Sf w ). Solving equations (7,8) and using the fact that the 
collective dynamics of the motors induces an active fric- 
tion [20] so that ico^ ac t (uj)5s u /2 = 8f u , we can obtain a 
linear relationship between Ssa, and F u and hence a fre- 
quency dependent correction to k e Q. The expression is 
complicated but simplifies in the experimentally relevant 
regime of stiff motors k m — > oo with the result 



1 



1 



z W L/ 2 0(l-0)r(G") 



fc effH fc cff(°) 



G> 



-3G'G" + 2(G') 



l\3 



(16) 



p-4 



Mechanical response of active gels 



where uj c ~ {CactG'L) . For w » w c we obtain an en- 
hancement of the effective stiffness of the elastic element: 



1 



1 



/c cff (oo) fc eff (0) 



1 



{G"f 



-3G"G" + 2{G'f . (17) 



G' 

v\3 



The effective low and high frequency elastic constants are 
shown in Fig. 3. Active cross- links always soften the zero 
frequency stiffness k e ff{0) of the elastic element, but at 
the same time always increase fc e //(oo) relative to the stiff- 
ness fco of a single filament. An estimate of the crossover 
frequency to c suggests that the stiffening may be relevant 
at the intermediate frequencies probed in experiments. 
From [20] we estimate ( act ~ ^p "^^ 1 ' , where N 



k Hff /k 




lVct /. 



Fig. 3: Plot of fc e ff(0) (red online) and fc e g-(oo) (blue on- 
line) for 1/km = 0, — = 0.4 as a function of L\Joo/k for 
Lpy/ao/n =0.2. 



is the cluster size, W\^ are the typical potentials, Wi,2 
typical motor conformation transition frequencies and I 
the motor's step length. From the expression for G above 
we estimate (at uq — 0) , G' = flf)t ,^ Tf ^ . Finally we obtain 



w, 



90^ 



90k B TL'i 
2 



^(^)(^) (l) 2 . Using ^4nm, 
VF 2 - Wi ~ 10fc s T, L p ~ 20 x 10 3 nm, L ~ 4/xm, A ~ 100, 
and wi^ — 100Hz, we obtain oj c ~ 10 _3 i?2;. 

We found that both the contractility and the stiffening 
of the active element are proportional to the square of the 
stall force of the motor cluster. This point and its impli- 
cation for the macroscopic shear modulus of a crosslinked 
network merit some discussion. 

Shear modulus of an active gel. First we 

use standard methods to relate the shear modulus of a 
crosslinked network of noninteracting elastic elements to 
the stiffness k e g of each elastic element [5,6]. We de- 
scribe the active gel as a cross-linked scmiflcxiblc poly- 
mer gel of monomer density p, where each monomer is 
a sphere of diameter a <C L p . The mesh-size is then 
£ 2 = I /{pa) [4-6]. Under a uniform shear, a point x 
of the network is deformed according to x — > x + a(x), 



with 7y = i(<9iOj + djtti) the applied strain. An elas- 
tic segment of orientation u and end-to-end distance Lq 
undergoes a relative change SLq/Lq = jijUiUj. This de- 
formation will in turn induce a tension o~ s = fc c g-<5Lo 
in the elastic unit, where l/fc e ff is the longitudinal re- 
sponse function of the unit. The corresponding contribu- 
tion to the stress tensor of the gel is <Ty = £,~ 2 (a s UiUj), 
where (...) denotes an average over the filaments' orienta- 
tion. For an isotropic filament distribution [uiUjU^ui) = 
{l/lb)[8ij8ki +8ikSji +5u6jk]- Defining the shear modulus 
E of the gel via a^ = 2E^ij [5,6], for an incompressible 
gel we obtain E = j^£_~ 2 L k e g. 

The presence of molecular motors acting as active cross- 
links modifies both the mean end-to-end length Lq of a 
filament strand and the stiffness k e /f of each elastic ele- 
ment. As shown earlier, an active crosslink stiffens fc e // 
at high frequencies {u> > uj c ), corresponding to time scales 
shorter than the relaxation time of the motor clusters, 
but longer than the relaxation time of a single filament. 
This correction, given in Eq. (17), is proportional to the 
square of the stall force /, which in turn is linearly pro- 
portional to the ATP activity A/i (the chemical potential 
difference between ATP and its hydrolysis products), for 
small A/i. This will therefore yield an active stiffening of 
order (A/i) 2 of the shear modulus of the network. To esti- 
mate the effect of activity on Lq we note that the network 
contains a bulk density n paB3 of passive crosslinks and a 
bulk density n mot of active crosslinks. Among the latter 
only a fraction r = k on /{k on + k ff) are bound, where 
k on and fc // are the motors binding and unbinding rates, 
respectively. Assuming only pairwise crosslinks and no 
dangling ends, the mean strand length between crosslinks 
is Lq = 7}pa/{rn mot + n pass ). The rates k on and k ff de- 
pend on A/i and arc finite at chemical equilibrium when 
A/x = 0. In general r is expected to depend linearly on 
activity for small A/i, i.e., r{Ap) s=s ro+riAp. This imme- 
diately gives Lq ~ Lq ' + A/iLg ' , and will therefore yields 
active corrections to the shear modulus linear in A/i. 

Expanding for small A/i, the shear modulus of the ac- 
tive gel will have the form -E activo = S paMive + ApE^ + 
{Ap) 2 E^ + O ((A/i) 3 ) . Explicit expressions for the vari- 
ous contributions can be obtained for instance in the limit 
of high frequencies using the expression (17) obtained ear- 
lier for k e ff{oo). The first term, S passive , is the shear mod- 
ulus of a passive cross-linked semiflexible gel [5,6]. The 
term linear in A/i comes from the change in the num- 
ber of crosslinks due to ATP consumption. The term 
quadratic in A/i has contributions from the active forces 
(calculated in this paper) , as well as smaller contributions 
from quadratic corrections to r. It can dominate at inter- 
mediate activities where it yields stiffening of the gel at 
high frequencies. At small A/i the correction due to the 
variation of the binding and unbinding rates with activity 
will dominate. The sign of this correction is controlled by 
the sign of T\ and is difficult to assess due to two com- 
peting effects. The unbinding rate k ff is known exper- 
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imentally and theoretically to increase with Afj, [14,21]. 
The binding rate k on is also expected to increase with A/i 
because once a motor cluster is bound to one filament, its 
directed motion along the filament allows it to explore a 
larger region of phase space and facilitates the binding to 
a second filament. This can also yield stiffening of the 
gel at very low activity with an elastic modulus increasing 
linearly with A/i if the increase of the binding rate with 
A/i dominates the increase in the unbinding rate (L < 
above). Detailed experiments are needed to address this 
question. 

Experiments have found that the addition of active 
crosslinkers such as myosin II can increase the shear mod- 
ulus of the network of several order of magnitudes im- 
plying that a quantitative comparison requires going into 
the non-pcrturbative regime. When the motor induced 
tension exceeds the bare tension, i.e., / > <tq, the parts 
of the filaments under compression will buckle and their 
response will be governed solely by the bending rigidity, 
k. Under such conditions, a complete calculation becomes 
more difficult. In the "high" frequency regime, however, 
the compressed parts of the filaments contribute negligi- 
bly to the force balance and the modulus is entirely con- 
trolled by the tense portion. A straightforward calcula- 
tion then yields E ~ 2^ 2 (rn mot /n pa3S )(K 1 / 2 / 3/2 / fc B T ) + 
0( s £ , -jfji ). When m mot ~ n paB3 , the shear modulus of the 

active gel scales as (//cro) 3//2 , in agreement with observa- 
tions [14, 17]. We stress that all our results only apply if 
the density of bound motors is larger than a critical den- 
sity required for a network of tense filaments to percolate 
trough the gel. For smaller values of rn mot , one expects 
little effect from the active crosslinks. Conversely, if the 
average number motor clusters bound to each elastic unit 
exceeds one, our formulae also breaks down as the por- 
tion of filaments between two successive motor clusters do 
not experience a large tension. These arguments suggest 
that both motor-induced contractility and stiffening will 
occur only in a narrow range of density of bound motors, 
in qualitative agreement with experiments. 

To summarise, we have studied a simplified microscopic 
model of a cross-linked active gel and shown that the non- 
linear elasticity and collective dynamics of the motors play 
an important role in the macroscopic mechanical proper- 
ties of the gel. In particular we show that clastic nonlinear- 
ities can lead to a gel which is contractile and stiffened by 
active elements above a characteristic crossover frequency 
due to the collective dynamics of the motors. 
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